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Figure 

1-1 
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Shock  velocity  (Us,  kn^'sec)  versus  particle  velocity  (up, 
km/sec)  for  grariite.  Solid  lines  represent  linear 
regressions  for  the  three  regimes.  Dashed  lines  represent 
the  calculations  for  the  mineral  mixture  model  with  the 
modal  composition,  and  the  dotted  line  represent  the 
calculation  for  the  oxide  mixture  model  of  constituent 
oxides.  These  calculations  do  not  suggest  any  significant 
difference  between  the  present  granite  and  the  Westerly 
granite  [McQueen,  1967].  The  elastic  velocity  obtained 
for  shot  #806  (in  parentheses)  is  considered  unreliable 
because  of  timing  uncertainties.  HEL  =  Hugoniot  elastic 
limit 

1-2  Pressure-density  relation  for  Hugoniot  states  of  granite.  9 

The  heavy  solid  curve  corresponds  to  the  Hugoniot 
equation  of  state  shown  in  Figure  1.  Solid  symbols 
represent  the  Hugoniot  states,  and  pluses  arc  the  measured 
partial  release  states  (Table  3).  Light  solid  lines  in  the 
present  study  serve  only  to  connect  the  data  sets  and  are 
not  intended  to  display  the  actual  release  path. 

1-3  Interface  particle  velocity  profile  during  shock  loading  of  a  10 

weathered  granite  under  dry  and  wet  conditions. 

1-4  Pressure-volume  relations  for  quartz  release  paths.  Data  22 

from  Podurets  et  al.  [1976  ],  Chhabildas  and  Miller 
[1985],  and  Swegle  [19%  ].  St  =  stishovite,  MP  =  mixed 
phase,  Qz  =  quartz.  FQ  =  fuzed  quartz. 

1-5  Pressure- volume  relations  for  quartz  release  paths.  Dotted  24 

and  broken  curves  illustrate  experimental  data,  and  solid 
curves  illustrate  Hugoniot  curves  and  calculated  release 
paths.  Solid  curves  are  computed  based  on  mixed 
Hugoniots  of  Qz  and  St  regimes  with  mass  fractions  at 
0.1  intervals  for  St.  Below  critical  pressures  (Table  7), 
the  release  paths  are  modeled  by  straight  lines.  The  zero- 
pressure  volume  of  post-shock  quartz  shocked  above 
some  50  GPa  was  taken  to  be  0.44  cmVg. 

1-6  Comparison  of  pressure-volume  relations  for  measured  27 

granite  release  paths  with  the  model  calculations.  HP- 
high-pressure,  LP  =  low  pressure,  MP  =  mixed  phase, 
and  EXj  =  diaplectic  glass. 
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1  -7  Post-shock  temperatures  for  quartz  and  granite  plotted  as  a  31 
function  of  peak  shock  pressure.  The  heavy  solid  line 
represents  preferred  values  for  post  shock  temperatures  in 
quartz  between  0  and  100  GPa.  Post  shock  temperatures 
calculated  for  granite  in  the  mixed  phase  region  are  shown 
as  the  solid  line.  Also  shown  are  experimental 
determinations  for  SiCh  in  both  the  quartz  and  stishovite 
stability  fields.  The  data  sets  of  Raikes  and  Ahrens 
[1979]  and  Boslough  [1988]  are  from  measurements  of 
post-shock  radiation.  The  data  of  Chhabildas  and  Miller 
[1985]  represent  calculated  values  determined  from 
measured  wave  profiles  and  the  Mie-Gruneisen  equation. 

The  dashed  line  is  the  calculation  of  Wackerle  for  SiC)2. 

1-8  Shock  temperatures  in  SiC)2  as  a  function  of  shock  32 

pressure.  The  heavy  solid  line  represents  preferred  values 
of  shock  temperature  in  the  low  pressure,  mixed  phase, 
and  high  pressure  regimes  of  Si02.  The  steep  P-T  slope 
calculated  for  the  mixed  phase  region  in  this  study 
contrasts  with  the  shallow  slope  calculated  by  Wackerle 
[1962].  The  data  of  Lyzenga  et  al.  [1983]  were  obtained 
from  pyrometer  measurements.  The  present  results  agree 
well  with  the  extrapolated  trends  from  this  data  set. 

1- 9  Shock  temperatures  in  granite  as  a  function  of  shock  34 

pressure,  ^fcrred  values  of  the  shock  temperature  of 
granite  in  the  low  pressure  and  mixed  phase  regions  are 
shown  by  the  heavy  solid  line.  In  the  mixed  phase 
region,  granite  temperatures  are  similar  to  those  in  Si02. 

The  temperature  at  the  top  of  the  mixed  phase  region  is 
“700  K  lower  than  calculated  values  for  the  high  pressure 
phase  of  granite  assuming  Yo  =  1  and  Cv  =  1.19  J/gmK. 

2- 1  Shock  velocity  (Us,  km/sec)  versus  particle  velocity  (Up,  42 

km/sec)  for  muscovite.  Shock  velocity  at  Up=0  is  within 
the  range  of  acoustic  bulk  sound  velocity  for  muscovite 
[Vaughan,  1986  ]. 

2-2  Pressure-density  relations  for  Hugoniot  state  of  43 

muscovite.  Solid  circles  are  from  present  study,  broken 
curve  is  for  mineral  mixtures,  and  dotted  curve  assumes 
oxide  mixture  model  [Telegin,  1980].  Pluses  are  for  the 
static  compression  for  muscovite  [Bridgman,  1949].  The 
zero-pressure  densities  at  points  1 ,  2  and  3  correspond  to 
mixtures  of  orthoclase  +  AI2O3+H2O,  of  wadeite  + 
kyanite  +  AI2O3+H2O,  and  KAlSi308  (hollandite)  + 
AI2O3+H2O,  respectively. 
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2-3  Pressure-temperature  projection  of  phase  boundaries  and  49 

muscovite  Hugoniot.  Circles  represent  the  Hugoniot 
calculation  using  the  heat  capacity  at  298K  and  squares 
represent  those  using  the  value  at  1000  K.  Abbreviations; 
Hol=hollandite  KAlSisOg,  Cbr=Al203;  Or=orthoclase 
and  Liq=liquid.  Solid  curve  was  experimentally 
determined  by  Huang  and  Wyllie  [1973  ]. 

2-4  Pressure-density  relations  for  muscovite  showing  52 

adiabatic  release  states.  Solid  circles  are  Hugoniot  states 
and  pluses  represent  measured  partial  release  states  (Table 
3).  Solid  lines  serve  only  to  connect  the  data  sets  and  are 
not  intended  to  represent  the  actual  release  path.  Arrows 
indicate  calculated  slopes  obtainol  from  eq.  (31). 
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OF  GRANITE 
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RI  02912 


1 


INTRODUCTION 

The  present  report  is  divided  into  two  chapters.  Chapter  1  deals  with  the  equation 
of  state  of  granite.  The  analysis  of  the  data  for  granite  utilizes  the  equation  of  state  of  mica 
(muscovite),  which  is  reported  in  Chapter  2.  Shock  Hugoniol  data  on  silicate  rocks  and 
minerals  all  demonstrate  major  shock-induced  phase  transformations  (e.g.  see  van  Thiel 
[1977];  Marsh  [1980];  Trunin  [1986]).  The  shock  response  of  this  class  of  materials  is  of 
importance  in  describing  shock-propagation  from  impact  and  explosive  sources  on  the  earth 
and  other  planetary  surfaces  as  well  having  application  in  the  study  of  the  Earth’s  interior. 

Granitic  rocks  are  the  most  common  component  of  the  upper  continental  crust,  and 
consist  mainly  of  quartz  and  feldspars.  The  Hugoniots  of  tectosilicates  and  rocks 
composed  chiefly  of  tectosilicates  are  all  quite  similar.  A  number  of  studies  have  been 
performed  to  measure  the  high-pressure  properties  of  these  minerals  and  their  mixtures 
(e.g.  Wackerle  [1962];  Ahrens  et  al.  [1969];  Grady  et  al.[1975];  Lyzenga  et  al.  [1983]). 
The  derived  equation  of  state  parameters  have  been  compared  with  theoretical 
considerations  (e.g.  McQueen  et  al.  [1963];  Ahrens  et  al.  [1969]).  The  onset  of  shock- 
induced  phase  transformations  usually  does  not  occur  at  the  thermodynamic  equilibrium 
pressure  or  temperatuic.  For  example,  in  the  case  of  quartz,  the  shock-induced  high 
pressure  phase  is  inferred  to  be  stishovite  [McQueen  et  al.,  1963]  and  this  identification  is 
consistent  with  shock  wave  recovery  experiments  (e.g.  Milton  and  DeCarli  [1963]).  For 
feldspars,  the  measured  Hugoniots  above  30  GPa  can  be  interpreted  in  terms  of  a  high- 
pressure  polymorph  with  the  hollandite  structure  [Ahrens  et  al.,  1969;  Sekine  and  Ahrens, 
1991]. 

Shock  recovered  samples  of  quartz  and  feldspars  from  above  25-30  GPa  indicate 
transformation  to  diaplectic  glasses  which  are  charaaerized  by  higher  refractive  indices  and 
densities  than  normal  fused  glasses  of  the  same  composition  [DeCarli  and  Jamieson,  1959; 
Heyman  and  Hdrz,  1990;  Kleeman,  1971;  Velde  et  al.,  1989;  Wackerle,  1962].  These 
diaplectic  glasses  have  been  intensively  studied  (e.g.  Stoffler  and  Homemann  [1972]).  Are 


these  diaplectic  glasses  characteristic  of  tectosilicates?  Are  they  produced  during  shock 
compression  or  during  isentropic  release?  If  formation  of  glass  occurs  upon  isentropic 
release,  what  is  the  phase  present  in  the  compressed  state  and  at  what  pressure  does  the 
shocked  tectosilicate  transform  to  diaplectic  glass? 

Recentiy,  isothermal  compression  of  tectosilicates  as  well  as  fayalite  in  the  diamond 
anvil  cell  at  room  temperature  indicates  that  transformation  to  a  dense,  amorphous  phase 
occurs  with  an  increase  in  the  coordination  of  silicon  and  aluminum  by  oxygen  [Hemley  et 
al.,  1988;  Williams  and  Jeanloz,  1989;  Williams  et  al.,  1990]. 

The  purpose  of  this  study  is  to  present  new  experimental  data  for  Hugoniot  and 
partially  released  states  of  a  weathered  granite,  to  synthesize  these  and  previous  data,  and  to 
generate  a  complete  equation  of  state.  Initially,  we  summarize  available  data  on  the 
equation  of  state  of  granite  as  well  as  its  major  constituents,  quartz  and  feldspar.  The 
release  adiabat  states  are  described  in  terms  of  the  third-order  Birch  Mumaghan  equation  of 
state  and  compared  with  those  for  quartz  and  feldspars.  We  then  construct  analytic 
expressions  for  the  release  isentropes  in  the  low-pressure,  mixed-phase,  and  high-pressure 
regimes.  The  mixed  phase  and  high-pressure  release  isentropes  are  described  in  terms  of 
frozen  release  isentropes  [Grady  et  al.,  1974;  Swegle,  1990]  and  upon  release  to  pressures 
in  the  2  to  10  GPa  range,  the  high-pressure  phase  is  assumed  to  transform  to  glass  as 
indicated  by  a  number  of  release  data  [Ahrens  and  Rosenberg,  1968;  Chhabildas  and 
Grady,  1984;  Chhabildas  and  Miller,  1985;  Podurets  et  al.,  1976], 

Experimental 

The  chemical  composition  of  the  weathered  granite  used  in  the  present  study  was 
determined  by  electron  microprobe  analyses  of  a  glass  formed  from  rock  powder  quenched 
in  water.  The  weathered  granite  was  melted  at  1600°C  in  air  to  obtain  a  homogeneous 
glass.  Moral  analysis  of  the  weathered  granite  was  determined  by  point  counting.  Table  2 
lists  the  chemical  composition  and  the  mode,  together  with  the  compositions  of  the  other 


3 


granites  for  which  equations  of  state  have  been  measured  [Chhabildas  and  Grady,  1984; 
Marsh,  1980;  Van  Thiel,  1977].  The  chemical  composition  of  all  the  granites  is  similar. 

The  present  shock  specimens  were  cut  from  a  disk,  10  cm  in  diameter  and  3  cm 
thick.  Bulk  densities  ranged  between  2.619-2.642  g/cm^,  while  crystal  densities  ranged 
from  2.640-2.645  g/cm^,  indicating  that  porosity  of  the  samples  is  low.  The  measured 
longitudinal  and  shear  sound  velocities  were  5.36  ±  0. 10  and  3.3  ±  0.2  km/sec, 
respectively . 

Shock  compression  was  conducted  by  launching  impactors  with  the  Caltech  25-mm 
bore  two-stage  light  gas  gun  and  a  40-mm  bore  propellant  gun.  Metal  flyer  plate-bearing 
projectiles  were  used  to  impact  the  samples.  The  metal  flyer  plates  were  A1  or  Ta.  In  each 
experiment,  the  impact  velocity  of  the  projectile  is  measured  by  the  flash  x-ray  method  and 
the  shock  wave  velocity  is  determined  by  measuring  the  travel  time  of  the  shock  wave 
through  the  sample  of  known  thickness  by  means  of  a  rotating-mirror  or  an  image- 
converter  streak  camera.  Rat  and  inclined  minors  were  employed  to  observe  the  shock 
wave  travel  time  [Ahrens,  1987].  The  Hugoniot  state  is  calculated  by  applying  the 
impedance  match  conditions  to  the  measured  initial  density,  impact  velocity  and  shock 
wave  velocity. 

A  partial  or  fully  released  state  is  also  determined  by  a  free-surface  or  buffer 
impedance  mismatch  observation,  reduced  by  way  of  the  Riemann-integral  formalism 
[Lyzenga  et  al.,  1983;  Rice  et  al.,  1958].  Lexan  and  polyethylene  were  used  as  buffer 
materials,  mounted  at  the  rear  of  the  sample.  Streak  camera  cut  offs  could  not  always  be 
observed  for  the  buffer. 

A  velocity  sensitive  interferometer  (VISAR)  [Baiker  and  Hollenbach,  1972]  was 
employed  to  investigate  compressive  wave  profiles  in  dry  and  water-saturated  granite  at 
low  stresses  (~10  GPa).  This  technique  makes  use  of  the  Doppler  shift  of  reflected  laser 
light  produced  by  the  motion  of  a  diffuse  surface.  Interference  fringes  proportional  to  the 
velocity  of  the  reflector  are  developed  in  a  modified  Michelsen  interferometer  and  recorded 
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using  photo-multiplier  tubes  (Burle  7764)  and  digitizing  oscilloscopes  (HP  541 1  ID).  The 
VISAR  used  in  the  present  experiments  incorporates  the  push-pull  modification  [Hemsing. 

1979]  for  improved  signal  quality.  The  time  resolution  of  our  VISAR  is  estimated  to  be  5 
ns  and  the  velocity  resolution  limit  is  about  10  m/s.  The  target  assembly  for  these 
experiments  consisted  of  an  aluminum  driver  plate  (1.5  mm  thick),  a  granite  target  (6.25 
mm),  a  diffuse  reflecting  aluminum  buffer  (0.75  mm),  and  a  LiF  window  (8  mm).  For  one 
shot  (#806),  the  target  was  held  in  an  aluminum  cup  and  was  water- saturated  in  vacuum. 
The  shock  equations  of  state  for  Al,  Ta,  LiF,  lexan  and  polystyrene  employed  in  these  and 
the  Hugoniot  experiments  are  summarized  in  Table  2.  The  travel  time  through  the  sample 
and  reflector  was  read  from  the  velocity  profile  and  corrected  to  calculate  the  velocities  of 
elastic  and  plastic  waves. 

RESULTS 

The  results  of  the  present  experiments  are  given  in  Table  3,  and  are  shown  in  Figs. 

1  and  2.  For  comparison,  the  Westerly  granite  data  [McQueen  et  al.,  1967]  are  also  plotted 
in  the  figures  and  included  in  the  analysis. 

The  shock-velocity  (Us)  versus  particle  velocity  (Up)  relations  shown  in  Fig.  1  can 
be  characterized  by  two  linear  regressions  (Us  =  Q)  +  SUp)  Co  =  4.70  km/sec  and  S  = 
0.46  in  a  range  of  Up  =  0.7  to  2. 1  km/sec  and  Cq  =  2.66  km/sec  and  S  =  1 .49  in  a  range  of 
Up  =  1.8  to  4.8  km/sec.  Moreover,  the  elastic  portion  of  the  Hugoniot  can  be 
approximated  by  a  linear  relation  of  Us  =  5.36  +  0.88  Up  (Up  =  0-0.6  km/sec).  The 
effect  of  water  saturation  on  elastic  shock  velocity  is  not  resolved  by  the  present 
experiments.  The  bulk  sound  velocity  resulting  from  the  measurements  of  Vp  and  Vs  for 
dry  granite  at  1  atm  is  3.77  ±  0.55  km/sec  which  is  close  to  Westerly  granite  data  [Marsh, 

1980] .  The  value  of  Vp  of  granite  increases  under  H20'Saturated  conditions  [Nur  and 
Simmons,  1969].  It  would  seem  that  the  presence  of  water  in  granite  has  a  critical 
influence  only  in  the  relatively  low  pressure  regime  (P  ^  0. 1  GPa).  It  appears  that  the 
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small  amount  of  water  which  can  be  injected  in  the  granite’s  pores  has  little  effect  in  the 
range  of  the  present  experiments. 

Table  4  summarizes  the  Us-Up  relations  for  various  granites.  Their  compositions 
are  given  in  Table  1.  Despite  their  chemical  variations,  the  Us-Up  relations  for  granites  are 
all  similar.  Qj  ranges  between  2.10  and  2.66  km/sec  and  S  is  between  1.49  and  1.63, 
when  Up  ranges  from  1.8  to  6  km/sec.  Above  Up  =  5  km/sec,  only  a  few  data  are 
available  (biotite-chlorite  granite  and  Soviet  granite)  and  Tmnin  et  al.[1988]  obtained  a 
linear  relation  for  the  Soviet  granite  (Table  4). 

Telegin  et  al.  [1980]  have  demonstrated  good  agreement  between  the  observed 
Hugoniot  and  the  calculated  Hugoniot  based  upon  the  oxide  mixture  model.  According  to 
this  model. 


C„  =  a„-(-a,p„-»-£  a^Zj 

i 

(la) 

S  =  b,+  b,p.+  X  b* 

(lb) 

where  Po  is  the  initial  density,  Zi  is  the  weight  percentage  of  component  oxide  i,  and  ao, 
ai,  ai,  bo,  bi,  and  bi  are  constants.  The  constants  reported  by  Telegin  et  al.  [1980]  were 
used  in  the  present  calculations.  The  calculations  are  restricted  to  the  high  pressure  phase 
region.  The  results  shown  in  Fig  1  (dotted  curve)  (Us  =  2.369  -i- 1.59  Up  km/sec) 
demonstrate  a  good  approximation  to  the  Hugoniot  data  in  the  high  pressure  regime. 

We  also  apply  a  mineral  mixture  model  [Al'tshuler  and  Sharipdzhanov,  1971]  to 
obtain  a  theoretical  Hugoniot  for  the  present  weathered  granite.  In  this  model. 


V(P)=Z  V(P)M_ 


(2) 


where  Vj  is  the  volume  of  constituent  mineral  i  at  pressure  P  and  Mj  is  a  mass  fraction  of 
mineral  i.  Using  the  Rankine-Hugoniot  equations.  Us  and  Up  are  computed  from  the 
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Shock  Velocity  (km/sec) 


0.0  1.0  2.0  3.0  4.0  5.0 

Particle  Velocity  (km/sec) 


Figure  1-1  Shock  velocity  (Us,  km/sec)  versus  particle  velocity  (up,  km/sec)  for 
granite.  Solid  lines  represent  linear  regressions  for  the  three  regimes. 
Dashed  lines  represent  the  calculations  for  the  mineral  mixture  model  with 
the  modal  composition,  and  the  dotted  line  represent  the  calculation  for  the 
oxide  mixture  model  of  constituent  oxides.  These  calculations  do  not 
suggest  any  significant  difference  between  the  present  granite  and  the 
Westerly  granite  [McQueen  et  al.,  1967].  The  elastic  velocity  obtained  for 
shot  #806  (in  parentheses)  is  considered  unreliable  because  of  timing 
uncertainties.  HEL  =  Hugoniot  elastic  limit. 
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resulting  P-V  relation.  The  granite  is  taken  to  be  a  mineral  mixture  consisting  of  the  modal 
composition  (Table  1).  The  Hugoniots  for  plagioclase  and  K-  feldspar  have  been 
recognized  to  be  insensitive  to  composition  and  the  microcline  Hugoniot  [Ahrens  et  al., 
1969;  Simakov  et  al.,  1974]  was  used  in  the  present  calculation.  The  quartz  Hugoniot  was 
taken  from  the  data  of  Wackerle  [  1962]  and  Trunin  et  al.  [  1970].  The  muscovite  Hugoniot, 
which  has  been  recently  determined  [Sekine  et  til.,  1991a],  was  also  employed  to  obtain  the 
calculation  of  the  weathered  granite  Hugoniot.  The  resulting  Us-Up  relation  is  shown  in 
Fig.  1  (dashed  line).  It  approximates  the  experimental  data  only  in  the  high-pressure 
regime  (Up  >  2.0  km/sec),  but  not  in  the  low  pressure  regime  (Up  <  2  km/sec). 

In  Fig.  2,  pressure-density  relations  are  plotted  for  the  present  granite  and  Westerly 
granite.  The  pressure-density  relation  may  essentially  be  divided  into  three  parts:  low- 
pressure  regime,  mixed  phase  regime  and  high-pressure  regime.  The  separation  between 
the  low-pressure  and  mixed  phase  regimes  is  not  complete  as  in  the  pressure-density 
Hugoniot  for  quartz  [Swegle,  1990]. 

The  goal  of  the  VISAR  experiments  was  to  determine  the  effect  of  water-saturation 
on  the  HEL  and  .shock  stmcture  in  the  double  shock-wave  regime.  Impact  conditions  are 
listed  in  Table  3  and  particle  velocity  profiles  for  the  aluminum-LiF  interface  are  shown  in 
Figure  3.  In  all  experiments,  two  wave  structure  is  visible,  but  it  is  less  well  defined  in  the 
case  of  the  wet  experiment.  The  shock  velocities  for  the  elastic  and  plastic  waves  were 
estimated  from  the  VISAR  profile  and  thickness  of  sample  by  assuming  that  the  elastic  and 
plastic  wave  velocities  of  Al-2024  are  constant  at  6.36  knVsec  [Marsh,  1980]  and  given  by 
the  Us-Up  relations  showi,  in  7  able  2,  respectively.  The  particle  velocity  in  the  Al  driver 
plate  is  approximated  by  taking  0.5  of  the  impact  velocity  and  that  in  tlie  reflector  Al  is  read 
from  the  VISAR  record.  The  calculated  results  are  lisicd  in  Table  3.  The  Us-Up  relations, 
which  are  included  in  Fig.  1 ,  show  good  agreement  with  the  Westerly  granite  data  but  a 
considerable  divergence  from  the  calculations  based  on  tJie  mineral  mixture  model.  This 
may  indicate  the  importance  of  material  strength  in  the  low  pressure  region. 


R 


Pressure  (GPa) 


2.5  3.0  3.5  4.0  4.5  5.0  5.5 

Density  (g/cm^) 


Figure  1-2  Pressure -density  relation  for  Hugoniot  states  of  granite.  The  heavy  solid 
curve  corresponds  to  the  Hugoniot  equation  of  state  shown  in  Figure  1. 
Solid  symbols  represent  the  Hugoniot  states,  and  pluses  are  the  measured 
partial  release  states  (Table  3).  Light  solid  lines  in  the  present  study  serve 
only  to  connect  the  data  sets  and  are  not  intended  to  display  the  actual 
release  path. 
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Velocity  (km/s) 


Time  (/^sec) 


Figure  1-3  Interface  particle  velocity  profile  during  shock  loading  of  a  weathered 
granite  under  dry  and  wet  conditions. 
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Table  1.  Chemical  Compositions  and  Modal  Compositions  of  Sevenil  Granites 
Used  for  the  Equation  of  State  Experiments 
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Table  2.  Equation  of  State  Parameters  of  Materials  Used  in  the  Present  Study 


Material 

Initial 

Density 

(g/cm3) 

Us= 

Co 

(krn/sec) 

=Co+S  Up  (km/sec) 

S 

Ref. 

A1  2024 

2.785 

5.328 

1.338 

1 

Ta 

16.656 

3.43 

1.19 

2 

LiF 

2.638 

5.15 

1.35 

1 

Lexan 

1.193 

2.449 

1.498 

1 

Polystyrene 

foam 

6.055 

0.243 

1.118 

1 

Ref 

1.  Marsh  [1980] 

2.  Mitchell  and  Nellis  [1981] 
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Table  4.  A  Comparison  of  Equations  of  State  for  Granites 
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The  HEL  (Hugoniot  elastic  lirmt)  has  been  determined  from  the  VISAR  records  and 
the  inclined  mirror  method.  The  latter  indicated  a  higher  value  of  Up  and  a  higher  pressure 
for  the  HEL  (Table  3).  This  may  be  due  to  the  fact  the  the  HEL  pressure  increases  slightly 
with  increasing  driving  shock  pressure  [Ahrens  et  al.,  1969].  The  HEL  wave  velocity 
demonstrates  some  scatter  for  the  three  shots  fired. 

In  Table  5,  the  final  shock  Hugoniot  states  of  granite  and  the  Al  reflector  are 
summarized  and  compared  with  the  measured  particle  velocity  of  the  reflector  at  the 
interface  of  the  LiF  window.  In  this  calculation,  the  Hugoniot  for  granite  was  assumed  to 
be  a  simple  linear  relation  (Us  =  3.77  +  1.52  Up)  connecting  the  two  points  of  the  lowest 
deformational  shock  state  and  the  bulk  sound  velocity  at  1  atm  in  the  Ug-Up  plot.  As  seen 
in  Table  5,  the  difference  in  the  calculated  and  measured  results  is  very  small  for  shots 
#805  and  #806.  For  shot  #803,  a  slight  optical  misalignment  resulted  in  non-normal 
incidence  of  the  laser  beam,  in  which  case  the  apparent  measured  velocity  is  less  than  the 
actual  particle  velocity. 

In  Figure  3,  both  the  elastic  and  plastic  portions  of  the  wavefront  are  characterized 
by  long  rise  times,  pjrobably  related  to  the  large  grain  size  of  the  samples.  Water-saturation 
results  in  a  further  increase  in  rise  time  and  obscures  the  inflection  point  between  the  elastic 
and  plastic  portions  of  the  wave.  According  to  Figure  3,  the  amplitude  of  the  plastic 
precursor  is  about  20%  lower  in  the  case  of  the  wet  shot.  The  lower  compressive  strength 
of  the  wet  sample  could  be  due  to  the  presence  of  water  along  cracks  reducing  the  effective 
normal  stress  and  therefore  the  frictional  strength  of  the  flaws,  thus  promoting  brittle 
failure. 

The  measured  release  states  of  the  present  granite  are  shown  in  Fig.  2,  together 
with  the  data  by  Swegle  [  1990].  The  release  of  shocked  granite  from  about  26-30  GPa 
Hugoniot  pressures  was  found  to  occur  along  pressure-density  paths  characteristic  of  the 
high-pressure  phase(s)  and  at  low  pressures,  final  unloading  states  were  close  to  the  low 
pressure  phase  initial  density.  Swegle  [1989]  also  observed  that  the  unloading  data  of 


15 


Table  5.  A  Comparison  of  VISAR  Data  With  the  Impedance  Match 
Calculations  for  Final  Shock  State 
Based  on  Equation  of  State  For  Granite 


Us=3.77+1.52  Up  km/s 


#805  Granite,  Hugoniot 
Reflected  from  A1 
Al,  Hugoniot 
Reflected  from  LiF 
Measured 
(Al-LiF  interface) 

#806  Granite,  Hugoniot 
Reflected  from  LiF 
Al,  Hugoniot 
Reflected  from  LiF 
Measured 
(Al-Lif  interface) 

#803  Granite,  Hugoniot 
Reflected  from  Al 
Al,  Hugoniot 
Reflected  from  LiF 
Measured 
(Al-LiF  interface) 


Us 

Up 

P 

(km/sec) 

(km/sec) 

(GPa) 

4.818 

0.689 

8.7 

4.958 

0.782 

10.2 

6.127 

0.597 

10.2 

6.098 

0.576 

9.8 

0  570 
(peak) 

5.002 

0.811 

10.6 

5.159 

0.914 

12.4 

6.274 

0.707 

12.4 

6.240 

0.682 

11.9 

.... 

0.661 

(peak) 

.... 

5.288 

0.998 

13.9 

5.464 

1.114 

16.0 

6.509 

0.883 

16.0 

6.469 

0.853 

15.4 

0.746 

(peak) 

Note:  The  upper  row  of  granite  corresponds  to  the  interface 
between  driver  plate  and  granite,  and  the  lower  one 
corresponds  to  that  between  granite  and  A 1  reflector.  The 
upper  row  of  A 1  corresponds  to  the  interface  between 
granite  and  Al  and  the  lower  respxands  to  that  between  Al 
’■'“flector  and  LiF  window. 
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granite  from  27  GPa  demonstrated  inirial  unloading  characteristic  of  the  high  pressure 
phase.  Swegle’s  [  1990]  data  indicate  that  the  reverse  transition  is  not  initiated  immediately 
upon  unloading.  This  behavior  can  be  seen  at  least  up  to  the  Hugoniot  pressure  of  54  GPa 
for  granite  [Swegle,  1990].  However,  the  release  path  observed  from  one  datum  at  the 
Hugoniot  pressure  of  124  GPa  indicates  a  shallower  slope  than  the  Hugoniot.  Swegle 
[1989]  provides  a  similar  unloading  profile  for  granite  at  92  GPa.  There  are  several  other 
data  which  indicate  that  the  release  path  slope  can  decrease  with  increasing  shock  pressure 
in  the  pressure -density  plot  (Swegle  [1990],  for  rhyolite). 


DISCUSSION 


The  shock  wave  data  are  used  to  obtain  equation  of  state  parameters  Kqs  (the  zero- 
pressure  adiabatic  bulk  modulus)  and  K'os  (its  first  pressure  derivative)  for  the  third  order 
Eulerian  finite  strain  (Birch-Mumaghan)  equation  of  state  for  granite.  We  use  the  approach 
of  Ahrens  and  Jeanloz  [1987];  Jeanloz  [1989];  Heinz  and  Jeanloz  [1984]  to  analyze  the 
Hugoniot  data. 

The  Eulerian  strain  /  is  expressed  as 


f  = 


1 

2 


-1 


(3) 


where  V  and  Vq  are  Hugoniot  and  ambient  volumes  of  the  sample  respectively.  The 
normalized  pressure  F hs  is  given  by 


(4) 


F„s=FHs  +  AF„  =  K„(l-2^.f3„+-) 

“  3f(l  +  2ff^[l  +  (2-1.5Y)f]*'” 

^  = _ PoTj^tr _ 

“  3f[l  +  (2-1.5Y)f] 

_f[l  +  (2-Y)f] 
l  +  (2-1.5Y)f 


(5) 

(6) 

(7) 


and 

^.  =  3{4-K:)/4  (8) 


where  y  is  the  Gruneisen  parameter  and  Ph  is  Hugoniot  pressure. 

The  parameters  Fhs.  fSH.  and  AFtr  are  obtained  from  the  experimental  data,  and 
the  parameter  Kos  and  K'os  are  obtained  from  a  linear  fit  to  equation  4.  y  is  calculated 
according  to  the  following  equations 


(9) 


and 


(10) 


where  a  is  the  volume  coefficient  of  thermal  expansion  and  Cp  is  the  specific  heat  at 
constant  pressure. 

Errors  in  Fhs  and  /3h  are  estimated  by  standard  error  propagation  techniques 
using  the  expressions  in  Ahrens  and  Jeanloz  [1987]. 


Low  Pressure  Rggimg 

For  granite  in  the  low  pressure  regime  (P  <  20  GPa),  Etr  =  0,  po  =  2.63  g/cm^,  Yo 
=  1.0  ±  0.5,  and  q  =  1.0  ±  1.0  were  used  in  the  present  computation.  These  values  cover 
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reasonable  ranges  for  Yo  and  q,  and  the  resulting  Kos  and  K'os  are  insensitive  to  the  values 
of  Yo  and  q  used.  An  unweighted  linear  least  squares  fit  to  the  experimental  data  below  a 
shock  pressure  of  20  GPa  including  the  Westerly  data  yields  Kos  =  57  GPa  and  K'os  = 

1.8.  A  Zero-order  fit  (K'os  =4)  of  the  experimental  data  yields  Kos  =  40  GPa.  Ultrasonic 
determinations  of  Kos  give  37  GPa  which  is  close  to  Kos  of  quartz. 

The  bulk  moduli  for  end-member  feldspars  range  between  67  and  106  GPa  [Angel 
et  al.,  1988].  The  value  of  Kos  =  57  GPa  for  granite  is  consistent  with  a  mixture  of  quartz 
and  feldspars.  The  equations  of  state  are  summarized  in  Table  6. 

High  Pressure  Regime 

Equation  of  state  parameters  are  calculated  for  granite  in  the  high-pressure  regime 
(P  >  39  GPa)  in  the  same  manner  as  for  the  low  pressure  regime.  The  values  of  Yo  £tnd  q 
are  not  well  constrained,  however.  Etr  is  estimated  to  be  0.8  ±  0.5  kJ/g  based  on  the 
transition  energies  of  quartz  to  stishovite  and  feldspars  to  the  hollandite  structures. 
Although  the  transition  energy  varies  with  feldspar  composition,  the  stated  error  range  for 
Etr  covers  the  whole  range  of  variation  of  feldspar  composition.  The  zero-pressure  density 
of  the  granite  in  the  high  pressure  regime  is  calculated  to  be  4.02  g/cm^  from  stishovite  and 
hollandite  densities.  The  hollandite  density  varies  also  with  the  initial  feldspar  composition 
and  was  approximated  by  the  modal  composition.  We  assumed  that  Yo  =  UO  ±  0.5  and  q 
=  1.0  ±  1.0.  With  these  parameters,  we  obtained  an  adiabatic  bulk  modulus  of  Kqs  =  304 
GPa  and  its  first  pressure  derivative  K'os  =  U2.  These  values  are  similar  to  estimates  for 
stishovite  based  on  quartz  Hugoniot  data  [Tan  and  Ahrens,  1990],  although  the  zero- 
pressure  density  differs  slightly  (by  5%).  The  datum  of  shot  #222  indicated  a  larger 
normalized  Fhs  with  4  times  greater  error  than  the  others.  Swegle  [  1 990]  gave  the  same 
equation  of  state  parameters  for  quartz  and  granite  in  both  the  low-pressure  and  high- 
pressure  regimes.  They  are  compared  in  Table  6.  Recent  analyses  of  Hugoniot  data  for 
various  feldspars  provide  ranges  of  Kos  K'os  as  well  as  the  zero-pressure  densities. 
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Table  6.  Summary  of  Equation  of  State  Parameters  of  Granite  (Gra),  Quartz  (Qz), 
and  Feldspar  (Fid)  for  the  I-X>w  Pressure  Regime  (LP)  and  High-Pressure  Regime  (HP) 
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based  on  shocked  feldspars  possessing  the  hollandite  structure  [Sekine  and  Ahrens,  1991]. 
The  equation  of  state  parameters  for  granite  are  between  those  of  quartz  and  feldspars  for 
both  the  low  and  high  pressure  regimes. 


Release  States 

The  Hugoniots  and  equation  of  state  parameters  for  quartz,  feldspars  and  granite  are 
quite  similar  as  summarized  in  Table  6.  The  release  adiabat  states  measured  for  shocked 
quartz  [Chhabildas  and  Miller,  1985;  Podurets  et  al.,  1976],  feldspars  [Ahrens  et  al.,  1969; 
Grady  and  Murri,  1976]  and  granite  [Sekine  et  al.,  1991b;  Swegle,  1989;  Swegle,  1990] 
also  show  similar  release  paths,  depending  mainly  upon  the  initial  shock  state.  These 
release  adiabats  have  been  measured  by  independent  methods:  impedance  mismatch  buffer 
method  [Ahrens  et  al.,  1969;  Podurets  et  al.,  1976;  Sekine  et  al.,  1991b],  manganin 
transducer  record  [Grady  and  Murri,  1976]  and  velocity  interferometer  (VISAR) 
[Chhabildas  and  Miller,  1985;  Swegle,  1989]. 

Quartz  has  been  studied  widely  under  dynamic  conditions.  Figure  4  shows  the 
release  adiabats  measured  for  quartz.  The  release  adiabats  from  shocked  quartz  in  the 
mixed  phase  regime  are  approximated  by  unloading  paths  in  which  the  quantity  of  low  and 
high  pressure  phase  in  the  mixture  is  frozen  down  to  about  8  GPa  [Grady  et  al.,  1974; 
Swegle,  1989].  At  this  point,  the  volumes  expand  to  the  initial  volume  with  further 
decreasing  pressure. 

The  release  adiabats  from  the  stishovite  regime,  however,  seem  to  occur  along 
paths  leading  to  less  dense  zero-pressure  material  than  the  zero-pressure  density  of  the 
high-pressure  form  (stishovite),  especially  when  the  release  pressure  reaches  within  the 
mixed  phase  regime.  It  appears  that  a  material  with  density  of  about  3.7  g/cm^  is  produced 
during  the  unloading  process  and  may  then  revert  to  a  diaplectic  glass  on  further  pressure 
release  to  ambient  pressure  and  temperatiue  as  observed  in  the  shock-recovery  experiments 
on  quartz.  This  apparent  density,  in  fact,  may  be  a  diaplectic  glass  with  six-fold 
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Figure  1  -4  Pressure-volume  relations  for  quartz  release  paths.  Data  from  Podurets  et 
al.  [1976],  Chhabildas  and  Miller  [1985],  and  Swegle  [1990].  St  = 
stishovite,  MP  =  mixed  phase,  Qz  =  quartz.  FQ  =  fuzed  quartz. 
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coordinated  silicon  with  oxygen  .  We  set  the  post-shock  quartz  density  to  be  2.27  g/cm^ 
[Stoffler  and  Homemann,  1972],  It  is  interesting  to  note  that  the  density  ratio  of  this 
proposed  material  to  stishovite  is  close  to  that  of  the  fused  silica  glass  to  quanz.  In  this 
context,  Chhabildas  and  Miller  [1985]  suggest  that  this  material  might  represent  a  high 
density  quartz  liquid  produced  above  the  melting  point  during  unloading. 

Figure  5  displays  calculations  of  release  adiabats  in  the  mixed  phase  regime  based 
on  the  mixing  of  two  Hugoniots  for  the  quartz  regime  and  the  stishovite  regime  in  a  range 
of  pressure  for  the  mixed  phase  regime.  In  this  model,  the  mass  fraction  of  stishovite  is 
frozen  in  the  mixed  phase  during  the  release  state  until  the  pressure  drops  down  to  a  critical 
value  represented  by  the  relation  P  (GPa)  =  32.2  -  85.5  V  (cm^/g)  shown  in  Fig.  5. 

Further  release  causes  the  formation  of  diaplectic  glass  and  the  final  volume  of  post-shock 
quartz  is  given  by  the  relation  V(cm3/g)  =  0.063M  +  0.377  at  ambient  conditions,  where  M 
is  the  mass  fraction  of  stishovite  at  the  Hugoniot  state,  as  listed  in  Table  7.  This  model 
calculation  is  compared  with  the  experimental  data  [Chhabildas  and  Miller,  1985;  Podurets 
et  al.,  1976;  Swegle,  1990].  For  release  paths  from  the  stishovite  Hugoniot  regime,  the 
zero-pressure  density  of  3.7  g/cm^,  Yq  =  1.0  and  q  =  1.0  were  taken  to  obtain  the  equation 
of  state  parameters.  The  release  data  points  were  taken  from  the  two  Hugoniot  pressures  of 
108  and  71  GPa  [Chhabildas  and  Miller,  1985].  The  results  are  listed  in  Table  7. 

As  seen  in  Fig.  2,  the  release  path  of  granite  from  30  GPa  also  appears  to  behave  in 
a  similar  manner  as  that  proposed  for  quartz.  Figure  6  compares  the  calculated  release 
paths  of  granite  in  the  mixed  phase  regime  with  the  experimental  data  The  calculations  are 
based  on  the  mixing  of  the  two  Hugoniots  for  the  low  pressure  and  high  pressure  regimes 
and  the  frozen  mass  fraction  of  the  Hugoniots  state  down  to  a  critical  pressure.  Below  the 
critical  pressure  given  by  P(GPa)  =  19.3  -  77  V  (cm^/g)  the  release  paths  are  represented 
by  straight  lines  on  which  the  formation  of  diaplectic  glasses  proceeds. 

By  applying  the  finite  strain  formalism  to  parts  of  the  release  adiabats  from  the 
Hugoniot  pressures  within  the  mixed  phase  regime,  we  obtain  the  K^s  and  K'qs  for  granite. 


Pressure  (GPa) 


Vo(DG) 


Figure  1-5  Pressure-volume  relations  for  quartz  release  paths.  Dotted  and  broken 
curves  illustrate  experimental  data,  and  solid  curves  illustrate  Hugoniot 
curves  and  calculated  release  paths.  Solid  curves  are  computed  based  on 
mixed  Hugoniots  of  Qz  and  St  regimes  with  mass  fractions  at  0.1  intervals 
for  St.  Below  critical  pressures  (Table  7),  the  release  paths  are  modeled  by 
straight  lines.  The  zero-pressure  volume  of  post-shock  quartz  shocked 
above  some  50  GPa  was  t^cn  to  be  0.44  cmVg. 
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Table  7.  Summary  of  Equation  of  State  Parameters  of  Model  Release  Adiabats  for  Quartz 
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The  results  are  listed  in  Table  8,  together  with  the  zero-pressure  density  of  granite,  the 
critical  pressure  and  density  by  release  path,  and  the  zero-pressure  density  of  post-shock 
granite.  We  assumed  the  same  density  of  2.27  g/m^  as  in  quartz  for  the  proposed  diaplectic 
glass  of  granite  [Stoffler  and  Homemann,  1972].  The  unloading  paths  below  the  critical 
pressure  are  associated  with  relatively  large  volume  changes  (Fig.  6).  The  critical  pressure 
(and  density)  may  vary  with  the  critical  pressure  of  the  Hugoniot  state  [Chhabildas  and 
Miller,  1985].  The  parameters  are  listed  in  Table  8.  For  the  release  path  from  the  high- 
pressure  regime,  the  equation  of  state  parameters  were  calculated  using  the  data  from  54 
GPa  Hugoniot  pressure  (Swegle,  1989].  The  results  are  listed  also  in  Table  8,  as  well  as 
the  zero-pressure  density  and  post-shock  granite  density  used. 

According  to  Swegle  [1989]  the  unloading  processes  of  granite  differ  from  those  of 
quartz.  In  crystalline  quartz,  the  reverse  transformation  is  well  described  by  the 
equilibrium  reverse  transition  and  by  frozen  phases  until  the  equilibrium  phase  boundary  is 
reached,  and  then  the  material  reverts  to  the  initial  phase.  This  behavior  is  also  the  case  for 
fine-grained  polycrystalline  quartz  [Grady  et  al.,  1974].  In  the  unloading  paths  of  granite, 
no  part  is  described  by  the  equilibrium  reverse  transition  although  they  are  characterized  by 
frozen  high  and  low  pressure  phases.  The  unloading  paths  of  feldspars  [Ahrens  et  al., 
1969;  Grady  and  Murri,  1976]  illustrate  a  different  behavior  from  that  of  quartz,  showing 
a  deviation  from  the  frozen  high  and  low  pressure  phase  mixtures  toward  the  low  density, 
well  above  the  expected  equilibrium  boundary  for  the  high  pressure  phase  (hollandite). 

The  conditions  for  formation  of  diaplectic  glasses  appear  to  be  restricted  to  only  part 
of  the  high  pressure  phase  present  in  the  Hugoniot  state  which  reverts  a  disordered 
diaplectic  glass  via  solid-solid  transition.  Shock-melted  glasses  have  been  also  recognized 
(c.g.  StSffler  and  Homemann  [1972]).  The  formation  of  shock-melted  glasses, 
characterized  by  closer  relation  of  the  physical  properties  to  normal  glasses  quenched  from 
melt,  requires  much  higher  shock  pressures  and  consequently  higher  shock  temperatures. 
The  release  adiabat  for  the  shock-melt  state  must  be  quite  different  from  the  release  adiabats 
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Figure  1  -6  Comparison  of  pressure- volume  relations  for  measured  granite  release  paths 

with  the  model  calculations.  HP-high-pressure,  LP  =  low  pressure,  MP  = 
mixed  phase,  and  EXj  =  diaplectic  glass. 
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discussed  here.  There  appears  no  significant  change  in  the  unloading  behavior  for  quartz 
release  paths  observed  up  to  the  Hugoniot  pressure  of  108  GPa  [Chhabildas  and  Miller, 
1985].  This  is  consistent  with  the  shock  temperature  measurements  of  quartz  to  support 
the  presence  of  crystalline  stishovite  without  melting  up  to  a  pressure  of  about  1 10  GPa 
[Lyzenga  et  al.,  1983]  although  the  shock  equation  of  state  for  quartz  at  much  higher 
pressures  reveals  some  phase  transformation  [Trunin  et  al.,  1970]. 

For  anorthite  glass,  the  adiabatic  releases  from  Hugoniot  pressures  up  to  40  GPa 
are  consistent  with  the  frozen  high-pressure  and  low-pressure  phase  mixture  [Boslough  et 
al.,  1986].  The  releases  from  much  higher  shock  pressure  than  90  GPa,  however,  showed 
a  significant  change  of  the  release  adiabat  slope.  It  seems  to  be  due  to  melting  according  to 
a  proposed  phase  diagram  for  anorthite  [Schmitt  and  Ahrens,  1983].  For  the  other 
feldspars,  no  significant  change  of  the  adiabatic  release  behavior  was  observed  at  least  up 
to  the  Hugoniot  pressure  of  46  GPa  [Ahrens  et  al.,  1969;  Grady  and  Murri,  1976].  The 
results  of  investigations  of  post-shock  samples  reveals  slightly  lower  pressure  values  for 
formation  of  shock-fused  glasses,  i.e.  about  43  GPa  for  feldspars  and  50  GPa  or  40  GPa 
for  quartz  [Stoffler  and  Homemann,  1972;  Tattevin  et  al.,  1990]. 

g-Y-E  JRglatiQns  for  Quaru  and  Granite 

The  pressure-volume  unloading  paths  depicted  in  Figures  4-6  provided  additional 
information  regarding  energies  and  hence  temperatures  achieved  during  tlie  shock  and 
subsequent  release  process.  This  offers  a  means  to  further  test  the  inferred  release  paths  by 
comparing  shock  and  post- shock  temperatures  with  available  data  and  theoretical 
expectations.  The  specific  energy  due  to  compression  is  obtained  from  conservation  of 
energy  across  the  shock  front: 

E„  =  P„(V-V)/2  (11) 

The  Hugoniot  energies  for  quartz  and  granite  are  listed  in  Tables  8  and  9  for  points 
in  the  mixed  phase  region  corresponding  to  increasing  fractions  of  frozen  high  pressure 


material.  Integration  of  the  P-V  release  paths  in  Figures  4-6  yields  the  total  energy  under 
the  release  isentrope,  Ef.  These  are  also  listed  in  Tables  8  and  9  for  quartz  and  granite, 
respectively.  Residual  (post-shock)  temperatures  were  calculated  from; 

X  =  \  +  (Eh  -E,)/Q  (12) 

where  Cy  is  the  average  specific  heat  at  constant  volume  [Robie  et  al.,  1978]  over  the 
temperature  interval  TrTo. 

The  resulting  post-shock  temperatures  for  quartz  and  granite  are  contained  in  Tables 
7  and  8  and  Figure  7.  For  quartz,  the  maximum  temperature  achieved  is  1673  K  which  is 
less  than  the  1  atm.  melting  temperature  of  Si02  (1996  K).  The  maximum  temperature 
achieved  in  granite  is  similar  (17 13  K  at  40  GPa).  We  therefore  infer  that  at  pressures 
between  15-40  GPa  both  quartz  and  granite  remain  solid  upon  isentropic  unloading.  As 
discussed  above,  Chhabildas  and  Miller  [1985]  infer  that  release  from  pressures  above  50 
GPa  results  in  a  high-density  quartz  liquid.  The  transition  energy  associated  with  the 
transformation  of  the  high  pressure  component  to  diaplectic  glass  can  be  estimated  using 
the  quartz  to  fused  quartz  transition  energy  (123  J/gm)  for  the  transforming  high-pressure 
phase  fraction.  The  effect  of  including  this  term  is  to  produce  a  slight  decrease  in  the 
calculated  post-shock  (and  shock)  temperature  (cf:  Figure  8). 

Also  shown  in  Figure  7  are  experimentally  determined  post-shock  temperatures  for 
quartz  [Raikes  and  Ahrens,  1979]  and  stishovite  [Boslough,  1988;  Chhabildas  and  Miller, 
1985].  The  present  results  are  in  good  agreement  with  the  highest  pressure  quartz  data  but 
are  -500  K  colder  than  the  stishovite  data  extrapolated  to  the  top  of  the  mixed  phase  region 
(-40  GPa).  Tlie  results  of  Chhabildas  and  Miller  [1985],  and  possibly  Boslough  [1988]  as 
well,  do  not  represent  complete  unloading  and  may  therefore  overestimate  the  residual 
temperature.  Tlie  calculation  of  Wackerle  [1962]  based  on  the  Mie-Gruneisen  Equation 
predicts  lower  temperatures  in  the  mixed  phase  region  and  a  PT  slope  in  the  stishovite  field 
which  is  much  steeper  than  suggested  by  the  experiments. 
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Postshock  Temperature  (K) 


Figure  1-7  Post-shock  temperatures  for  quartz  and  granite  plotted  as  a  function  of  peak 
shock  pressure.  The  heavy  solid  line  represents  preferred  values  for  post 
shock  temperatures  in  quartz  between  0  and  100  GPa.  Post  shock 
temperatures  calculated  for  granite  in  the  mixed  phase  region  are  shown  as 
the  solid  line.  Also  shown  are  experimental  determinations  for  SiC>2  in  both 
the  quartz  and  stishovite  stability  fields.  The  data  sets  of  Raikes  and  Ahrens 
11979]  and  Boslough  [1988]  are  from  measurements  of  post-shock 
radiation.  The  data  of  Chhabildas  and  Miller  [1985]  represent  calculated 
values  determined  from  measured  wave  profiles  and  the  Mie-Gruneisen 
equation.  The  dashed  line  is  the  calculation  of  Wackerle  for  Si02. 
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Shock  Temperature  (K) 


Pressure  (GPa) 


Figure  1-8  Shock  temperatures  in  Si02  as  a  function  of  shock  pressure.  The  heavy 
solid  line  represents  preferred  values  of  shock  temjjerature  in  the  low 
pressure,  mixed  phase,  and  high  pressure  regimes  of  Si02.  The  steep  P-T 
slope  calculated  for  the  mixed  phase  region  in  this  study  contrasts  with  the 
shdlow  slope  calculated  by  Wackerle  [1962].  The  data  of  Lyzenga  et  al. 
[1983]  were  obtained  from  pyrometer  measurements.  The  present  results 
agree  well  with  the  extrapolat^  trends  from  this  data  set. 
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The  temperature  decrease  along  the  isentropic  path  ftom  the  shock  state  to  complete 
unloading  can  be  used  to  infer  the  Hugoniot  temperamre.  Using  the  Mie-Gruneisen 
equation  of  state  together  with  the  assumption  that  yA^  is  constant  results  in  the  following: 


T„  =  T,exp 


7(v, -V„) 


(13) 


where  the  subscript  H  refers  to  Hugoniot  conditions,  o  to  the  initial  state,  and  r  to  the  post¬ 
shock  state.  Hugoniot  temperatures  calculated  in  this  manner  are  displayed  in  Tables  7  and 
8  and  Figures  8  and  9.  The  incorporation  of  the  quartz  -  bised  quartz  phase  transition 
energy  has  only  a  small  effect  (Figure  8).  Temperatures  in  the  mixed  phase  region  are 
found  to  increase  much  more  rapidly  with  pressure  than  those  in  either  the  low  pressure 
phase  or  high  pressure  phase.  This  is  due  to  large  increases  in  energy  caused  by  the  large 
volume  change  associated  with  the  transformation.  It  is  interesting  to  note  that  the 
calculated  temperature  for  Si02  at  the  top  of  the  mixed  phase  region  is  very  similar  to  the 
extrapolated  value  for  stishovite  at  that  pressure  based  on  the  measurements  of  Lyzenga  et 
al.  [1983].  This  contrasts  with  the  calculation  of  Wackerle  [1962]  who  predicted  a  very 
shallow  FT  slope  in  the  mixed-phase  region  but  a  much  steeper  slope  than  observed  in  the 
stishovite  region.  Figure  9  demonstrates  that  shock  temperatures  for  quartz  and  granite  are 
very  similar  in  the  mixed-phase  region.  The  calculated  temperanue  for  granite  at  the  top  of 
the  mixed-phase  region  is  ~7(X)  K  higher  than  that  calculated  for  the  high  pressure  phase  at 
nearly  the  same  pressure.  Temperatures  for  the  high  pressure  phase  were  obtained  using 
the  approach  of  McQueen  et  al.  [1970]  and  the  classical  limit  of  3R  for  the  specific  heat. 
These  results  demonstrate  that  the  present  model  for  both  quartz  and  granite  are  consistent 
from  an  energetic  viewpoint  and  produces  temperatures  in  reasonable  accord  with  expected 
behavior. 
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Shock  Temperature  (K) 


Pressure  (GPa) 


Figure  1-9  Shock  temperatures  in  granite  as  a  function  of  shock  pressure.  Preferred 
values  of  the  shock  temperature  of  granite  in  the  low  pressure  and  mixed 
phase  regions  are  shown  by  the  heavy  solid  line.  In  the  mixed  phase 
region,  granite  temperatures  are  similar  to  those  in  SiC)2-  The  temperature  at 
the  top  of  the  mixed  phase  region  is  -TOO  K  lower  than  calculated  values  for 
the  high  pressure  phase  of  granite  assuming  y©  =  1  and  Cv  =  1.19  J/gmK. 
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CONCLUSION 

New  equation  of  state  data  for  granite  have  been  obtained  and  combined  with  the 
Westerly  granite  data  of  McQueen  et  al.  [1967].  the  shock  equation  of  state  can  be 
described  by  two  linear  relations,  Us  =  4.70  +  0.46  Up  (km/sec)  for  a  range  of  Up  from 
0.7  to  2  km/sec  and  Us  =  2.66  +  1.49  Up  for  a  range  of  Up  from  2  to  5  km/sec.  These 
shock  equations  of  state  are  comparable  with  other  data  previously  determined  for  various 
granites.  The  third  order  Birch-Mumaghan  equation  of  state  parameters  are  Kos  =  57  GPa 
and  K’os  =1.8  for  the  low-pressure  regime  (po  =  2.63  g/cm^)  and  Kos  =  304  GPa  and 
K'os  =  1.2  for  the  high-pressure  regime  (po  =  4.02  g/cm^).  These  parameters  lie  between 
those  of  quartz  and  feldspar  in  the  low-pressure  regime,  and  between  those  of  stishovite 
and  the  hoUandite  form  in  the  high-pressure  regime. 

For  the  release  adiabats  from  the  mixed  phase  regime,  the  mixed  Hugoniot  model  is 
applied  and  compared  with  experimental  data.  In  this  model,  the  release  adiabats  are 
calculated  for  the  high-pressure  and  low-pressure  phase  mixture  with  the  same  mass 
fraction  as  at  the  Hugoniot  state  to  critical  pressures  given  by  P(GPa)  =  19.3  -  77  V 
(cm3/g). 

Further  releases  below  the  critical  pressures  are  associated  with  large  volume 
changes  related  to  formation  of  diaplectic  glasses  as  quenched  modifications  of  the  high- 
pressure  fractions. 

Based  on  the  measurements  of  partially  released  states  of  granite  as  well  as  quartz 
shocked  in  the  high-pressure  regime  (P  >  50  GPa),  it  is  proposed  that  the  high-pressure 
form  relaxes  isentropically  to  a  metastable,  inteimediate  phase  characterized  by  dense, 
highly-disordered  material.  This  material  is  subsequently  quenched  to  diaplectic  glass. 

The  process  of  diaplectic  glass  formation  differs  from  that  of  shock-fused  glass  which  can 
result  only  from  much  higher  Hugoniot  pressure.  The  release  adiabat  associated  with  the 
shock-fused  glass  formation  must  be  distinguished  from  the  release  adiabat  discussed  in  the 
present  study.  The  inferred  release  curves  from  the  mixed  region  for  quartz  yield  Hugoniot 
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temperatures  which  are  considerably  higher  than  those  of  Wackerle  [1962],  but  are 
consistent  with  the  data  of  Lyzenga  et  al.  [1983].  Moreover,  the  post  shock  temperatures 
predicted  agree  closely  with  the  data  of  Raikes  and  Ahrens  [1979]  at  low  pressures  and  the 
data  of  Chhabildas  and  Miller  [1985],  and  Boslough  [1988]  at  higher  pressures.  These 
agreements  lend  further  support  to  the  release  model  developed  here. 
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INTRODUCTION 

Hugoniot  equation  of  state  of  minerals  and  rocks  provide  the  basis  of  describing  shock 
wave  propagation  ftom  intense  explosions  in  the  Earth,  the  eflPects  of  meteorite  impact  on 
the  Earth  and  the  planets,  as  well  as  contributing  to  our  knowledge  of  the  Earth's  interior 
and  the  accretion  of  the  Earth.  Small  amounts  of  water  in  the  interior  can  be  requestered  in 
hydrous  minerals  such  as  serpentine,  brucite,  amphiboles  and  micas.  These  hydrous 
minerals  may  play  a  critical  role  in  the  evolution  of  the  interiw,  as  well  as  in  the  physical 
and  chemical  nature  of  the  mantle.  Muscovite  is  also  among  the  potassium-bearing 
minerals  that  could  influence  the  chemistry  of  mantle. 

Since  there  are  no  previous  equation  of  state  data  available  for  the  most  common  crustal 
hydrous  mineral,  muscovite,  it  was  the  subject  of  the  present  study.  The  present  data 
complement  recent  other  shock  wave  studies  of  the  hydrous  minerals  [Duffy  et  al.,  1991; 
Tybuiczy  et  al.,  1990]. 

Previously,  Cummings  [1968],  Short  [1968b],  Short  [1968a],  Horz  and  Ahrens  [1969], 
and  Lambert  and  Mackinnon  [1984]  investigated  the  pwst-shock  response  of  muscovite 
from  labwatory  experiments  and  from  matwial  in  the  vicinity  of  underground  explosions. 
They  observed  characteristic  features  such  as  kinking,  mosaic  extinction,  melting,  and 
vesculation.  These  field  and  laboratory  studies  lacked  detailed  knowledge  of  the  equation 
of  state  of  muscovite.  Moreover,  studies  of  the  equation  of  state  of  granite  rocks  and  their 
syntheses  require  shock  wave  date  for  constituent  minerals. 

EXPERIMENTAL 

The  natural  muscovite  of  the  present  study  was  from  Methuen  Township,  Ontario 
(Harvard  University  Mineralogical  Museum,  Cat  #1 12791).  The  composition  is  listed  in 
Table  1 .  The  apparent  size  is  about  40x50x2mm.  The  average  crystal  density,  determined 
by  Archimedean  method  was  2.834610.0026  g/cm^.  This  is  in  good  agreement  with  x-ray 
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Table  1.  Chemical  Composition  (wt  %)  of  muscovite 


h. 

£ 

d 

£ 

Si02 

45.87 

45.20 

45.41 

45.26 

47.30 

AI2O3 

38.69 

38.46 

38.52 

38.40 

36.05 

Na20 

0.64 

0.60 

0.62 

“ 

0.21 

K2O 

10.08 

10.44 

10.91 

11.82 

10.77 

H2O 

4.67 

4.64 

4.53 

4.52 

4.9 

Total 

99.95 

99.34 

99.99 

100.00 

99.23 

a.  Hurlbut,  C.  S.  Jr.  [1956]  Remainder  -  not  listed  is  0.10  wt 

%MgO. 

b.  Eugsteretal.[1972]  Remainder  -  not  listed  are  FeO, 

FeiOs,  MnO,  Li203,  Ti02,  etc.  and  summed  to  be  about 
0.4  wt.% 

c.  Calculated  from  present  idealized  muscovite  (Ko.92Nao.08) 

Al3Si3Oi0(OH)2. 

d.  Calculated  from  pure  moscovite  KAl3Si3Oi0(OH)2. 

e.  Vaughan  and  Guggenheim  [1986]. 
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density  measured  for  a  muscovite  with  a  similar  composition  [Vaughan  and  Guggenheim, 
1986].  The  entire  sample  was  imbedded  in  epoxy  resin  and  cut  into  rectangular  aliquots  (7 
X  13  X  2mm). 

Shock  loading  was  conducted  with  the  Caltech  25mm  two-stage  light  gas  gun  using 
metal  flyer  plate  bearing  projectiles  to  impact  samples  at  speeds  of  up  to  6.3  km/sec. 
Projectile  velocity  just  prior  to  impact  was  determined  by  double  exposure  flash  x-ray 
photography  using  two  15  nsec  flash  x-ray  sources  and  electronic  time  interval  counters 
[Jeanloz  and  Ahrens,  1977],  Shock  wave  velocities  in  the  target  and  buffer  materials  were 
determined  by  recording  the  destruction  time  of  the  mirrors  through  an  image  converter 
streak  camera  and  xenon  light  source  [Ahrens,  1987].  A  40-mm  propellant  [Ahrens  et  al., 
1969]  was  used  to  accelerate  lexan  projectiles  fitted  with  either  Cu  or  Ta  flyer  plates  to 
velocities  up  to  2.4  km/sec.  Projectile  velocities  were  measured  using  flash  x-ray  images. 
Sample  and  buffer  mirror  shock  velocities  were  obtained  from  observations  of  shock- 
induced  loss  of  mirror  reflectivity  at  successive  reflective  surfaces  as  recorded  by  a 
rotating-mirror  streak  camera. 

The  flyer  and  driver  plates  used  in  this  study  were  the  same  to  obtain  the  symmetric 
impact  conditions.  The  assumed  equation  of  state  parameters  are  listed  in  Table  2.  The 
particle  velocity  behind  the  shock  front  and  pressure-density  states  were  calculated  through 
the  impedance  match  method  and  the  Rankine-Hugoniot  equations,  respectively. 
Uncertainties  were  determined  by  standard  error  propagation  techniques  [lackson  and 
Ahrens,  1979]. 

Partial  releasee  states  were  obtained  by  measuring  the  shock  wave  transit  time  through 
low  impedance  buffers  in  contact  with  the  sample.  The  pressure  and  particle  velocity  at  the 
buffer-sample  interface  were  determined  from  the  measured  shock  velocity  and  the  known 
equation  of  state  for  the  buffer  material  (lexan).  An  upper  bound  for  the  density  of  the 
partial  release  state  can  be  calculated  by  integrating  the  Riemann  integral  over  a  linear 
pressure-volume  (P-V)  path  [Lyzenga  and  Ahrens,  1978].  For  all  of  the  shots,  streak 
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camera  cutoffs  could  be  observed  for  the  lexan  arrival  mirrors  mounted  on  the  rear  sample 
surface. 

In  several  experiments,  polystyrene  foam  buffers,  as  well  as  lexan  mirrors,  were 
mounted  on  the  rear  sample  surface  to  observe  partial  release  states.  Table  2  summarized 
the  shock  wave  equation  of  state  parameters  for  flyers,  drivers,  and  buffers  used  in  the 
present  study. 

RESULTS  AND  DISCUSSION 

Present  experimental  results  and  calculations  are  listed  in  Table  3  for  the  Hugoniot  and 
partial  release  states  determined. 

Hugoniot  States 

Peak  shock  pressures  achieved  ranges  between  20  and  142  GPa.  The  relationship 
between  the  shock  velocity  (Us)  and  particle  velocity  (Up)  is  presented  in  Fig.  1.  The  data 
set  can  be  described  by  a  linear  regression  of  the  seven  points: 

Us=Co  +  S  Up  (7) 

where  the  Co  is  fitted  to  be  4.63  ±0.12  km/sec  and  the  constant  S  is  to  be  1.27  ±0.04.  The 
Co  determined  is  the  same  within  experimental  errors  as  the  average  value  measured  by 
Brillouin  scattering  technique  for  a  single  crystal  muscovite,  of  which  composition  is 
similar  to  ours  [Vaughan  and  Guggenheim,  1986].  The  parameters  Co  and  S  do  not 
change  significantly  (C©  =  4.62  ±0.12  and  S  =  1.27  ±  0.04)  upon  also  employing  the 
acoustic  value  of  Co  =  4.56  ±  0.40  to  the  data  set.  The  zero-pressure  bulk  modulus  (Kos) 
and  the  first  pressure  derivative  (Ko's)  can  be  calculated  using: 


^S“PoCo^ 

(8) 

^o's=4S-l 

(9) 

They  are  60.5±3.1  GPa  and  4.1±0.2,  respectively. 

Pressure-density  states  attained  in  the  shock  compression  of  muscovite  are  depicted  in 
Fig.  2,  and  compared  with  the  model  calculations  based  upon  the  assumed  Hugoniots  for 
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Figure  2-1.  Shock  velocity  (Us,  km/sec)  versus  particle  velocity  (Up,  km/sec)  for 
muscovite.  Shock  velocity  at  Up=0  is  within  the  range  of  acoustic  bulk 
sound  velocity  for  muscovite  [Vaughan  and  Guggenheim,  1986]. 
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Pressure  (GPa) 


Density  (g/cm®) 


Figure  2-2.  Pressure-density  relations  for  Hugoniot  state  of  muscovite.  Solid  circles 
are  from  present  study,  broken  curve  is  for  mineral  mixtures,  and  dotted 
curve  assumes  oxide  mixture  model  [Telegin  et  al.,  1980].  Pluses  are  for 
the  static  compression  for  muscovite  [Bridgman,  1949].  TTie  zero-pressure 
densities  at  points  1,  2  and  3  correspond  to  mixtures  of  orthoclase  + 
AI2O3+H2O,  of  wadeite  +  kyanite  +  AI2O3+H2O,  and  KAlSisOg 
(hollandite)  AI2O3+H2O,  respectively. 
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Table  2.  Shock  equation  of  state  parameters  for  flyers,  drivers,  and  buffers  used  n  the  present  study 


Material 

Initial  density 
(g/cm3) 

Us-Up  relation 
(km/sec) 

Measured  range 
of  Up  (km/sec) 

Ref. 

Cu 

8.92 

Us=3.91+1.51  Up 

0.2~4.2 

1 

Al,  1100 

2.712 

Us=5.38+1.34  Up 

1.5-5. 1 

1,2,3 

Ta 

16.656 

Us=3.43+1.19  Up 

0.8-4.9 

1 

lexan 

1.193 

Us=2.419+1.321  Up 

3.2-5.2 

1 

Us=2.449+ 1.498  Up 

0.4-28.0 

polystyrene  foam 

0.055 

Us=0.243+1.118  Up 

0.2-5 

1 

1.  Marsh  [1980] 

2.  Mitchell  and  Nellis  [1982] 

3.  Brown  and  Shatter  [1984] 
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non-reactive  mineral  mixture  (e.g.  KAlSi30g  (orthoclase  to  30  GPa  and  hollandite  above 
30GPa)  +Al203-^H20)  [Al’tshuler  and  Sharipdzhanov,  1971]  and  oxide  mixture  (e.g. 
K20-HAl203-HSi02+H20)  [Telegin  et  al.,  1980].  In  the  mineral  mixture  model 
calculations,  the  total  volume  (V)  at  a  given  pressure  is  calculated  to  be 

V(P)  =  ImiVi(P)  (10) 

i 

where  mi  is  a  fraction  of  the  individual  constituent  mineral  i  with  a  volume  of  v  at  the 
pressure  P.  The  Hugoniots  for  KAISi308,  AI2O3,  and  H2O  can  be  found  in  the  literature 
[Ahrens  et  al.,  1969;  Marsh,  1980;  Mitchell  and  Nellis,  1982;  Simakov  et  al.,  1974], 

Since  the  Hugoniots  for  alkali  feldspars  and  plagioclase  are  insensitive  to  composition,  the 
corresponding  feldspar  to  the  present  muscovite  can  be  taken  to  be  microcUne  for  which 
Hugoniot  has  been  determined  up  to  191  GPa  [Simakov  et  al.,  1974].  In  Figure  2,  the 
static  compression  dates  [Bridgman,  1949]  are  added  to  describe  the  low  pressure  region. 
However,  the  shock  wave  data  appear  to  indicate  greater  compressibility  than  the  static 
data.  Formal  fits  for  equation  of  state  parameters  from  the  4  GPa  static  data  yields  the 
values  of  Kos=  108+5  GPa  and  Ko's=22±4. 

Telegin's  model  of  oxide  mixture  Hugoniot  can  be  applied  to  synthesize  the  muscovite 
data  in  an  approximate  range  of  particle  velocity  from  2  to  5  km/sec  [Telegin  et  al.,  1980]. 
This  corresponds  to  a  pressure  range  of  40  to  about  150  GPa.  For  multi-component 
systems. 


Us=Co(po,  Z)  S(Po,  Z)  Up  (11) 

where  parameters  Co(po.  Z)  and  S(po,  Z)  are  given  from  the  initial  density  (p^)  of  rocks 
and  weight  percentage  (Z)  of  oxide  constituents  in  rocks. 
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(12) 


Co=  ao  +  aiPo+IZjai 

i 

S=  bg  +  b{Po  +  ]LZj[bj 


(13) 


where  constants  a©,  ai  ai  and  bo,  bi,  bi  have  been  obtained  statistically  [Telegin  et  al.. 


1980]. 

As  seen  in  Fig.  2,  the  measured  densities  are  in  good  agreement  with  the  two  model 
calculations  in  the  pressure  range  of  about  80  to  140  GPa. 

The  shock  wave  data  are  reduced  to  an  isentrope  by  ^plying  the  third  order  Birch- 
Mumaghan  equation  of  state.  A  formalism  is  the  finite  strain  model  to  analyze  the 
Hugoniot  data  [Ahrens  and  Jeanloz,  1987;  Heinz  and  Jeanloz,  1984;  Jeanloz,  1989]. 
The  Eulerian  strain  f  is  expressed  as 


f 


1^ 

2 


fp.'“ 


(14) 


where  p  and  po  are  Hugoniot  and  ambient  densities,  respectively.  The  normalized  pressure 
Fhs  is  given  by 


l-^[(l  +  2f)^^^-l 

3f(l  +  2f)^^^[l  +  (2-1.5y)f] 


Ph 


(15) 


where  Y  is  the  Gruneisen  parameter  and  Ph  is  the  Hugoniot  pressure. 


fo  ^9 


Y  =  YcJ 


where  Yo  = 


aK, 


os 


Cp  p( 


(16) 


a  is  the  volume  coefficient  of  thermal  expansion  and  Cp  is  the  specific  heat  at  constant 
pressure.  The  value  of  a  for  muscovite  is  unknown  and  we  use  a=35xl0‘^  K'l  from 
muscovite  [Guggenheim  et  al.,  1986].  Cp  is  assumed  to  be  Cv=0.818  J/g  k  which  is  given 
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by  Robie  et  al.  [1978].  By  using  Kos=60.5±3.1  GPa  from  the  shock  wave  equations  of 
state,  the  value  of  0.72±0.04  is  calculated  for  Yo-  The  q  values  of  1.011.0  were  taken  in 
the  present  analyses.  Without  a  phase  transformation  the  internal  energy  of  transformation, 
Etr  =  0,  and  a  linear  least  squares  fit  using  Eq.  (17)  to  the  combined  data  set,  yields  the 
parameters  of  Kos=52±14  GPa  and  Ko’s=3.2±0.3. 


Ff,s  =  FHs  +  AF„  =  K„(l-2y3H+...)  (17) 


where 


^  t:  \ 


V  '’oj 


3f[l  +  (2-1.5Y)f] 


(18) 


and 


f3H  = 


f|l-^P-Y)f| 

l  +  (2-1.5Y)f 


(19) 


The  present  shock  wave  data  of  muscovite  do  not  indicate  any  strong  evidence  for  a 
rapid  transformation  to  a  high-pressure  phase  in  the  Us-Up  or  pressure-density  plots. 
However,  the  present  data  cannot  exclude  the  phase  transitions  which  are  associated  with 
small  volume  changes.  There  are  several  possible  decomposition  reactions  for  muscovite 


(KAl3Si3Oi0(OH)2)  as  follows: 

KAl3Si3Oi0(OH)2=KAlSi3O8(Sa)+Al2O3+H2O  (20) 

KAl3Si30 1 0(OH)2=KAlSi3O8  •H20(SaHy)-t-Al203  (2 1 ) 

KAl3Si3Ol0(OH)2=l/2K2Si4O9(Wad)-i-Al2SiO5(Ky) 

+  l/2Al203-tH20  (22) 

KAl3Si3Oi0(OH)2=KAlSi3O8(HoI)+Al2O3+H2O  (23) 


where  Sa=sanidine,  SaHy=sanidine  hydrate,  Wad=wadeite  [Kinomura  et  al.,  1975], 


Ky=kyanite,  and  Hol=ho!landite  (Ringwood  et  al.,  1967).  The  calculated  volume  changes 
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associated  with  the  above  reactions  are  1 1.92,  -0.27,  -1 1.52,  and  -24.64  cm^/mol  at  the 
zero-pressure,  respectively.  The  initial  densities  are  shown  in  Fig.  2  for  the  decomposition 
products.  Reaction  (20)  is  limited  to  high  temperatures  and  relatively  low  pressures 
[Chatteijee  and  Johannes,  1974],  and  reaction  (21)  is  also  limited  to  relatively  low 
temperatures  and  pressures  [Seki  and  Kennedy,  1964]. 

To  understand  the  stability  of  muscovite,  a  phase  diagram  was  constructed  using 
available  thermodynamic  data  for  muscovite  and  its  decomposition  products. 
Thermodynamic  data  are  taken  from  Robie  et  al.  [1978]  and  are  extrapolated  outside  their 
applicability  to  muscovite  above  1000  K.  Thermodynamic  functions  for  H2O  were  taken 
from  Halbach  and  Chatteijee  [1982].  Gibbs'  free  energy  values  for  wadeite,  K2Si409, 
and  KAlSisOg  (hoUandite)  were  computed  for  the  following  equilibrium  reactions: 

2KAlSi308(Sa)=K2Si409(Wal)-«-Al2Si05(Ky)+Si02(Coes)  (24) 

and 

K2Si409  (Wad)+Al2Si05(Ky)+Si02(St)=2KAlSi308(Hol)  (25) 

We  approximated  the  boundary  conditions  as  follows  [Liu  and  Bassett,  1986]: 

P(GPa)=6.3+0.002T('’C)  for  reaction  (24)  (26) 

and 

P(GPa)=8.3-K).0035T  (°C)  for  reaction  (25)  (27) 

where  the  effects  of  both  temperature  and  pressure  on  volumes  were  neglected  as  a  first 
a[^roximation.  The  resultant  equilibrium  conditions  for  reaction  (23)  are  expressed  as 

P(GPa)=7.840.0047T  ("C)  (28) 

which  is  shown  in  Fig.  3. 

Temperatures  along  the  Hugoniot  were  calculated  using  the  method  of  McQueen  et  al. 
[1967].  We  integrated  the  equation: 


dV  IvJ  2Cv 


(29) 
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Pressure(GPa) 


Figure  2-3.  Pressure-temperature  projection  of  phase  boundaries  and  muscovite 
Hugoniot.Circles  represent  the  Hugoniot  calculation  using  the  heat  capacity 
at  298K  and  squares  represent  those  using  the  value  at  1000  K. 
Abbreviations:  Hol=hollandite  KAlSi308,  Cor=Al203;  Or=orthoclase  and 
Liq=liquid.  Solid  curve  was  experimentally  determined  by  Huang  and 
Wyme[I973]. 
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and 


dS  ^(Vq- 

dV  2T 

The  resulting  shock  temperatures  are  shown  in  Figure  3.  The  muscovite  Hugoniot  is 
expected  to  intersect  the  boundary  defined  by  reaction  (23)  about  8  GPa.  However,  the 
present  Hugoniot  data  obtained  for  muscovite  indicate  that  muscovite  remains  present  at 
least  up  to  about  60  GPa,  where  the  Hugoniot  temperature  is  about  1000°C.  If  muscovite 
decomposes  at  a  lower  shock  pressure,  the  Hugoniot  is  expected  to  coincide  with  that  of 
the  mineral  mixture  model  above  30  GPa  because  the  constituent  KAlSisOg  transforms  to 
the  hollandite  structure  [Ahrens  and  Liu,  1973;  Ahrens  et  al.,  1969]. 

The  equilibrium  pressure  few  reaction  (22)  involving  wadeite  phase  has  a  lower  pressure 
than  reaction  (23).  Since  there  is  no  indication  that  transition  wadeite  occurs  in  the 
microcline  Hugoniot  data,  it  is  hard  to  discuss  the  reaction  (22)  in  the  intermediate  pressure 
range. 

The  equation  of  state  parameters  of  muscovite  (Kos=52±4  GPa,  Ko’s=3.2±3)  are 
comparable  with  those  of  brucite  [Duffy  et  al.,  1991]  (Kos=51±4  GPa,  Ko’s=5.0±0.4) 
and  serpentine  [Tyburezy  et  al.,  1990]  (Kos=63.5±6.7  GPa;  Ko  s  =2.75+0.62)  which 
have  been  recently  determined.  The  elasticity  measurement  of  muscovite  by  Brillouin 
scattering  technique  gave  the  average  bulk  sound  velocity  of  4.56±0.40  km/sec,  [Vaughan 
and  Guggenheim,  1986]  which  corresponds  to  a  bulk  modulus  of  59.0±5.2  GPa.  This 
value  is  consistent  with  our  results.  A  finite  strain  analysis  of  the  tremolite  Hugoniot  data 
below  about  65  GPa  [Simakov  and  Trunin,  1980]  yields  the  Kos=76±6  GPa  and 
Ko  s=5.1±0.5.  The  static  compression  data  of  portlandite  Ca(OH)2,  which  is 
isomorphous  with  brucite  yields  Kos=37.8±1.8  GPa  and  Kos -5.210.7  at  room 
temperature  [Meade  and  Jeanloz,  1990].  It  is  noteworthy  that  these  hydrous  minerals  vary 
greatly  in  their  water  content  from  about  2  wt  %  of  tremolite  to  about  30  wt  %  of  brucite, 
but  that  they  have  similar  equation  of  state  parameters. 


V)  +  P 

-  (30) 
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Partial  Release  States 


Measured  partial  release  states  for  muscovite  are  listed  in  Table  3  and  depicted  in  Fig.  4. 
Although  the  shock  equation  of  state  for  lexan  (buffer  material)  has  been  determined  only  in 
the  Up  range  up  to  about  5.2  km/sec  [Marsh,  1980],  we  used  a  single  linear  regression  of 
the  high  pressure  portion  to  calculate  the  observed  release  states.  The  release  path  observed 
from  the  lowest  Hugoniot  pressure  about  20  GPa  displays  an  irregular  shape.  Taking  into 
account  the  experimental  errors  (Table  3),  the  shock  wave  may  have  been  attenuated  in 
polystyrene,  causing  us  to  measure  too  low  a  shock  velocity  in  the  buffer.  The  partial 
release  states  obtained  from  the  Hugoniot  pressures  below  1 10  GPa  indicate  that  the 
adiabatic  release  paths  are  steeper  or  almost  have  the  same  slopes  in  the  Hugoniot  density- 
pressure  plane.  These  steep  paths  have  been  interpreted  as  indicating  retention  of  the 
denser  high-pressure  phase  during  the  release  (e.g.  Swegle  [1990]).  As  seen  in  Fig.  4, 
however,  the  release  paths  from  the  Hugoniot  pressures  of  130  to  140  GPa  indicate  slightly 
shallower  paths  than  those  from  the  lower  pressures.  The  shallower  paths  mean  the 
expected  devolatilization  during  release.  The  local  slqje  of  isentropic  expansion  of  a  point 
on  the  Hugoniot  can  be  approximated  [McQueen  et  al.,  1967] 


^dP'l  _  Po(2Kos-Ph7Ph) 
UpJh  Ph{7Ph-{2  +  Y)Po} 


(31) 


The  resulting  slope  obtained  shows  a  simple  increase  in  the  slope,  as  seen  in  Fig.  4.  The 
uncertainty  of  the  calculated  slope,  however,  becomes  greater  with  increasing  pressure  and 
it  may  be  hard  to  say  that  the  initial  slope  of  isentropic  expansion  is  steeper  than  that  of 
Hugoniot  at  high  pressures. 
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Pressure  (GPa) 


Density  (g/cm®) 


Figure  2-4.  Pressure-density  relations  for  muscovite  showing  adiabatic  release  states. 

Solid  circles  are  Hugoniot  states  and  pluses  represent  measured  partial 
release  states  (Table  3).  Solid  lines  serve  only  to  connect  the  data  sets  and 
are  not  intended  to  represent  the  actual  release  path.  Arrows  indicate 
calculated  slopes  obtained  from  eq.  (31). 
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According  to  the  investigaticHis  of  micas  in  experimentally  shocked  gneiss  [Lambert  and 
Mackinnon,  1984]  localized  melting  of  muscovite  begins  upon  release  fiom  some  30  GPa 
and  is  completed  at  70  GPa.  Note  that  the  shock-melting  was  congruent  and  the 
decomposition  products  of  muscovite  were  not  observed  experimentally.  Moreover,  the 
final  shock  states  in  the  shock-recovery  experiments  were  achieved  via  multiple 
reverberation  while  those  of  the  equation-of-state  experiments  were  determined  under  a 
single  shock  compression.  The  final  pressures  in  recovery  experiments  are  those  in  the 
stainless  steel  containers  and  the  initial  shock  pressure  can  be  calculated  as  low  as  10  to  15 
GPa. 

Phase  relations  of  muscovite  under  static  pressure  conditions  [Huang  and  Wyllie,  1973] 
indicate  that  musco-'ite  decomposes  into  orthoclase.  corundum,  and  vapor  below  about  1 
GPa  and  825  °C,  and  that  it  will  melt  incongruently  above  that  pressure  and  temperature.  If 
muscovite  melts  incongruently  upon  release  fiom  the  shock  compression,  then  the  structual 
water  in  the  muscovite  may  be  expelled  at  considerably  low  pressures  during  the  pressure 
release  process.  The  pressures  determined  as  partially  released  states  are  apparently  too 
high  to  result  in  volatilization.  When  the  pressure  of  partial  release  state  comes  down  to  the 
~1  GPa  pressure  level,  we  expect  fallow  release  paths  and  much  lower  densities.  This 
behavior  is  observed  in  the  case  of  brucite  [Duffy  et  al.,  1991]. 

Conclusions 

(1)  Shock  equation  of  state  data  of  muscovite  are  determined  between  20  and  140  GPa, 
and  can  be  fit  by  a  linear  Us-Up  relation: 

Us=4,62+1.27  Up  (km/sec). 

(2)  Muscovite  apparently  remains  stable  along  its  HugonioL  Third  order  Birch- 
Mumaghan  equation  of  state  parameters  are:  Kos=52±4  GPa  and  Kos  =3.2±0.3.  These 
parameters  arc  comparable  with  the  other  hydrous  minerals  such  as  brucite,  serpentine,  and 
trenwlite,  although  the  water  content  varies  from  2  to  30  wt  %  among  these  minerals. 
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parameters  are  comparable  with  the  other  hydrous  minerals  such  as  brucite,  serpentine,  and 
trcmolite,  although  the  water  content  varies  from  2  to  30  wt  %  among  these  minerals. 

(3)  Thermodynamic  and  theoretical  calculations,  however,  suggest  that  muscovite 
dehydrates  into  the  hollandite  KAlSi308,  AI2O3,  and  H2O  with  only  a  small  volume 
change  above  a  pressure  near  70  GPa. 

(4)  Observed  release  paths  indicate  a  change  of  the  slope  in  the  density-pressure  plane  near 
80  GPa,  suggesting  vaporization  upon  adiabatic  pressure  release. 
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